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Abstract 

It is shown that there exists a possibility of amplification of amplitudes of quan- 
tum pion modes with low momenta in the late time of chiral phase transition by using 
the Gaussian wave functional approximation in the 0(4) linear sigma model. It is also 
shown that the amplification occurs in the mechanism of the resonance by forced oscil- 
lation as well as the parametric resonance induced by the small oscillation of the chiral 
condensate. These mechanisms are investigated in both the case of spatially homoge- 
neous system and the spatially expanded system described by the Bjorken coordinate. 
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§1. Introduction 



One of the recent interests in the context of the relativistic heavy ion collisions is to 
clarify the nature of matter at very high energy density. Especially, it is interesting to 
investigate the dynamics of the chiral phase transition in connection with the problem of the 
formation of a disoriented chiral condensate (DCC). 1 ) -10 ) Recently, Ikezi, Asakawa and one 
of the present authors (Y.T.) showed that there was a possibility of the formation of DCC 
by taking account of both the quantum fluctuations around the chiral condensate and the 
mode-mode coupling of quantum meson modes. n ) In the formation of DCC, it is necessary 
that the low momentum components of the inhomogeneous chiral condensate were grown. 
Some authors pointed out that the parametric resonance was seen when the formation of 
DCC occurred. 12) ' 13) 

Similarly, the amplitudes of quantum meson modes with low momenta may be grown 
in the chiral phase transition as well as the low momentum components of chiral conden- 
sate. 14 )' 15 ) In order to investigate the time evolution of the chiral condensate as an order 
parameter of chiral phase transition and of the fluctuation modes around it, it is suitable to 
adopt the time-dependent variational approach to dynamics of quantum fields in the 0(4) 
linear sigma model in terms of a squeezed state or a Gaussian wave functional. 16 )' 17 ) This 
method is based on the time-dependent variational method with a Gaussian wave functional 
by using a functional Schrodinger picture. 18 ) 

In this paper, we will show that the quantum meson modes are amplified even in the 
late time of chiral phase transition in both the cases of uniform and spatially expanding 
system in the 0(4) linear sigma model. Then, it will be pointed out that the possible 
mechanism of amplification is the resonance mechanism by the forced oscillation as well as 
the parametric resonance which was seen in the DCC formation. The unstable region in 
the parameter space is also presented concretely. This paper is organized as follows. In the 
next section, we recapitulate the time-dependent variational approach with a Gaussian wave 
functional to the 0(4) linear sigma model. The time evolution of chiral condensate and the 
quantum meson modes are governed by the Klein-Gordon type equation of motion and the 
Liouville-von Neumann type equation of motion, respectively. The numerical results are also 
given in §3. In §4, the time evolution of quantum meson modes in the late time of chiral 
phase transition is investigated. In both the case of the uniform and the spatially expanding 
system, it is shown that there exists a possible solution of equation of motion for quantum 
meson modes, which reveals a parametric resonance and/or resonance by forced oscillation 
induced by the small oscillation of the chiral condensate. In §5, the unstable regions are 
depicted for quantum meson modes, in which the amplitudes of meson modes are amplified 
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by the mechanism of parametric resonance and/or the resonance of forced oscillation. The 
last section is devoted to the summary. 

§2. Recapitulation of Gaussian wave functional approach to 0(4) linear sigma 

model 

In this section, we derive the equations of motion for the chiral condensate and quantum 
meson modes in the 0(4) linear sigma model. We apply the time-dependent variational 
method based on the functional Schrodinger picture to the 0(4) linear sigma model. 

2.1. Time- dependent variational approach with a Gaussian wave functional 

Following Refs.16) and 17), we derive equations of motion for the chiral condensate and 
quantum fluctuations. We start with the following Hamiltonian : 

/I 1 m 2 A 

d*x{-nl(x) + -iy Va {x)f + — Va {xf + -{ Va {xff - h<po(x)} , (2-1) 

where a runs 0~3. The index indicates the sigma field and 1~3 indicate the pion fields. 

We adopt the following Gaussian wave functional as a trial wave functional in the frame- 
work of the functional Schrodinger picture : 

^[v?(£c)]=Arexp(i(7r|v9-^)-(^-^l^-^l^-^ , (2-2) 

where Af is a normalization factor and G ab (x,y,t), E ab (x,y,t), if a (x,t) and n a (x,t) define 
the real and imaginary parts of the kernel of the Gaussian and its average position and 
momentum. Here, we have used simple notations as 

(n\(p) = / d 3 x^2n a (x,t)ip a (x) , 

J a=0 

3 

(<p\G\(p)= f[d 3 xd 3 yJ2 Va (x)G ab (x,y,t) Vb (x) , (2-3) 

a,b=0 

and so on. In the functional Schrodinger picture, the operator n a (x), which is a conjugate 
operator of the field operator v? (aj), is regarded as —i8/8ip a (x). The expectation value (O) 
for the field operator O is easily calculated such as 

((fi a (x)) = <p a (x,t) , 
(n a (x)) = 7t a (x,t) , 

(ip a (x)ip b (y)) = (p a {x, t)(p b (y, t) + G ab (x, y, t) , 

(7i a (x)n b (y)) = n a {x,t)7r b (y,t) + ^G-^x^^) + AE ■ G ■ E ab {x,y,t) , (2-4) 
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where we have used the following shorthanded notation : 

3 

A-B ab (x,y,t) =X!/ d 3 zA ac (x,z,t)B cb ( y z,y,t) . 

c=(T 

It is understood that (p represents the mean field which is identical with the chiral condensate 
in this model, and G represents quantum fluctuations around the mean field. 

The time dependence of the Gaussian wave functional is determined by the time de- 
pendence of variational functions G, E, (p and 7f, which is governed by the time-dependent 
variational principal : 

dt ( i §i~ A ) = °- (2 ' 5) 

The variational equations give the canonical equations of motion for (<p a , n a ) and (G ab , S ab ) 
with the Hamiltonian (H). As a result, we obtain the following equations of motion : 

= 7T , 

tt = |(V 2 - m 2 - ±<p 2 - ^trG(x, x)) - ^G(x, a;) J y? + hS a0 , (2-6) 
G = 2(E-G + G-E) , 

t = - V 2 + 2S 2 + i(-V 2 + m 2 + A (v?2 + t rG)) + ±{(pCp + G) . (2-7) 

These are a set of basic equations of motion for the chiral condensate <p and quantum 
fluctuations G. 

2.2. Reformulation of equations of motion for quantum fluctuations 

In this subsection, we reformulate the equations of motion for quantum fluctuations in 
order to investigate clearly the time evolution of quantum meson modes. This reformulation 
has been firstly carried out in Ref.17). We introduce the reduced density matrix M. ab '■ 

M ab (x,y,t) + -8 (x-y) = ■/*f\~f\\]> ( 2 ' 8 > 

2 y (7i a (x)7i b (y)} i(Tr a (x)ip b (y)) J 

where (p a = <p a — <p a and n a = n a — n a . Thus, the reduced density matrix can be expressed 
in terms of G and E as 

, A ( -2iG-E G \ 

M = , , . 2-9 

ViG" 1 +AE-G-S 21E-G) 

The time dependence of this reduced density matrix is governed by the Liouville-von Neu- 
mann type equation : 19 ^ 

iM ab (x, y, t) = [H, MUix, y, t) , (2-10) 
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where the generalized Hamiltonian 7i ab has particularly simple form as 

S ab \ x3 . 



H ab (x,y,t)=\ " \5\x-y) (2-11) 

r ab (x,t) o 



with 



r ab (x,t) = -V 2 5 ab ± M 2 b (x,t) , (2-12) 
M 2 b (x, t)= (m 2 + ^f(x, t) ± ^trG(a;, x, t)jd ab 

+ ^<fi a (x,t)<f b (x,t) + ^G ab (x,X,t) . 

Equations of motion derived from (2-10) are equivalent to (2-7). 

From the structure of Ai, it can be checked that the reduced density matrix satisfies 
a relation M 2 = 1/4. The eigenvalues of the reduced density matrix are thus ±1/2. We 
introduce an eigenvector of A4 ab , (u a ,v a ), with eigenvalue 1/2, i.e., 

v b (y,t) J 2 yv a (x,t) J 

Then, («*, — v*) is also an eigenvector with eigenvalue —1/2. The eigenvectors are conve- 
niently normalized to 

ldMu { :\x,t)% { r\^t)+^\x,tru^\x,t)) = ±5 mn p-u) 

a=0^ 

with a ± sign for eigenvalues ±1/2. For the particular normalization condition we have 
adopted the eigenvectors u a and v a provide the following spectral decomposition for the 
reduced density matrix A4 ab - 

u { a\x,t) 



M ab (x,y,t) = ^[(^[^^(^(l/,*)* ^\y,ty) 



. U ( a\x,t) 



(2-15) 



where sum runs over eigenstates with eigenvalues ±1/2 only. From Eqs.(2-9) and (2-15), the 
Gaussian kernel G and £ can be expressed in terms of u and v, for example, 

G ab (x,y,t) = • (2T6) 

n 

Thus, u a can be regarded as mode functions for quantum meson modes. 
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The equations of motion for mode functions can also be written by using the generalized 
Hamiltonian matrix : 

id (uaM\( S ab \fu b (x,t)\ 

\v a (x,t) J \r ab ( x ,t) J \v b (x,t) J 

The normalization condition (2-14) is preserved by these equations. Using the explicit form 
of the mean field operator r ab , we see that the mode function u a satisfy the set of coupled 
Klein-Gordon-type equations 

(nS ab + M 2 b (x,t))u b (x,t) = . (2-18) 

The matrix elements M 2 b are found to be diagonal in the isospin index, i.e.,M^ 6 = M 2 5 ab . 
Let us assume that y?o 7^ and (p 1 = </? 2 = <^ 3 = 0, that is, the chiral condensate points in 
the sigma direction. Then the Eqs.(2-6) and (2-18) read 

{□ + M 2 (x,t) - ^ (x,t) 2 ^ <p (x,t) = h , (2-19) 
(n + M 2 a (x,t))u a (x,t) = . (2-20) 

It should be noted here that, if the condensate is uniform spatially with translation 
invariance, namely, the condensate depends on only time, it is possible to carry out the 
Fourier transformation for u a (x,t) as follows : 



u a (x, t) = -= J d 3 ku k a (t)e^ . (2-21) 

Then, the above forms are substituted into Eq.(2-20) and <p depends only on time, we found 
following equations : 

d 2 t + M 2 (x,t) - ^<fo(t) 2 ) <p (t) = h , (2-22) 



3 

{d? + k 2 + M 2 a (x,t)}u*(t)=0. (2-23) 

These equations are investigated later. 

2.3. Geometry for the spatially expansion 

In this subsection, we derive the equations of motion that is taken account of the geometry 
of one-dimensional spatial expansion in z-direction. 20 -* For this purpose, the convenient 
variables, namely, the proper time r and usual rapidity variable rj are defined by 

r = Vt 2 ~^, V=l^l^- (2-24) 
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for one dimensional scaling case. 21 -* D'Alembertian in Eqs.(2-19) and (2-20) is rewritten by 
using these variables : 

□ = d 2 T + -d T - \d 2 - d 2 ± , (2-25) 

where <9 2 = d 2 + d 2 . 

If the condensate does not depend on (x±,rj), namely, the condensate depends on only 
proper time r, the quantum fluctuation modes u a can be expressed as follows: 7 ^ 



u a (x±,v,r) = -= [u^(r)e ik ^ ( 



ik ± -xj_ 



(2-26) 



where we have defined J k = Jd 2 k ± dk v /(2n) 3 ^ 2 . Then, the equations of motion (2-19) and 
(2-20) are recast into the following equations : 



(d 2 + \d T ) + M 2 (r) - ^(r) 



(p (r) = h , 



d 2 T + 



K + 1/4 



+ k 2 ± + M 2 (x ± ,r],T) 



«M(r) = 0, 



where momentum k v appears in the combination {k 2 + l/4)/r 2 . 



§3. Numerical calculation 



(2-27) 



(2-28) 



Let us demonstrate qualitatively the time evolution of the mean field, which is in the 
sigma-direction only, and of quantum meson mode functions in the case of an uniform system 
in Eqs. (2-22) and (2-23) without spatially expansion. We assume <pi — with % — 1 ~ 3. 
In the numerical calculation, we adopt the box normalization with the spatial length being 
L in each direction. We then impose the periodic boundary conditions for the fluctuation 
modes, namely the allowed values of momenta are k x = (2-7T / L)n x and so on, where n x 
is integer. The fluctuation modes labeled by (n x ,n y ,n z ) are included each direction up to 
n 2 = nl+Uy+n 2 , = 8 2 . This corresponds to the momentum cutoff A ~ lGeV(990 MeV) since 
we have adopted the collisional region as L 3 = (10 fm) 3 . Further, we assume that fluctuation 
modes of the pion fields are identical one another, which are denoted as u\ = u\ = = v%. 

In Fig.l, the time evolution of chiral condensate is depicted. In order to avoid the 
complexity of problems with respect to the initial conditions in the relativistic heavy ion 
collisions, we only demonstrate the time evolution qualitatively with <p (t = 0) = 40 MeV 
and (po = 0. It is shown that the chiral condensate approaches to the vacuum value gradually 
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Fig. 1. The time evolution of the chiral con- 
densate is depicted. The horizontal and 
vertical axes represent time and the value 
of chiral condensation, respectively. 




Time{fm/c) 

Fig. 2. The time evolution of the fluctuation 
mode with n = in the cr-direction is de- 
picted. 




Fig. 3. The time evolution of the fluctuation 
mode with n = in the 7r-direction is de- 
picted. 




Fig. 4. The time evolution of the fluctuation 
mode with n = 1 in the 7r-direction is de- 
picted. 



Fig. 5. The time evolution of the fluctuation 
mode with n = 2 in the 7r-direction is de- 
picted. 



Fig. 6. The time evolution of the fluctuation 
mode with n = 1 in the cr-direction is de- 
picted. 



with oscillation. In Fig. 2, the time evolution of the amplitude of quantum sigma meson mode 
with n = 0, that is Uq =0 , is shown. This mode function only oscillates. On the other hand, 
as is seen in Fig. 3, the amplitude of the quantum pion mode with n = 0, that is u^ =0 , is 
amplified. The amplitude of pion mode with n — 1 is also amplified weakly as is seen in Fig. 4. 
However, as is seen, for example, in Fig. 5, the amplification is not realized in the modes 
with n > 2. Also, the sigma meson modes with n > 1 as well as n = are not amplified 



8 



as is seen in Fig. 6. Thus, we conclude that, in this parameterization, the amplitudes of the 
lowest and the first excited pion modes are amplified, but other quantum meson modes are 
not grown. 

Next, let us demonstrate qualitatively the time evolution of mean field and of quantum 
meson mode functions in Eqs.(2-27) and (2-28) with spatially expansion. We assume the same 
conditions with respect to the chiral condensate as those in the case of the homogeneous 
condensate, namely <po 7^ and (pi — for % — 1 ~ 3, and the same initial condition for r 
instead of t. We then impose the periodic boundary conditions for the fluctuation modes, 
k x = (2ir/L)n x , k y = {2^/L)n y and k n = (27r/a)n v , where n x ,n y ,n v are integer, L is spatial 
length in x-y direction and a is dimensionless parameter. The fluctuation modes labeled 
by (n x , n y , n n ) are counted up to {^f) 2 {n 2 x + n 2 y ) + (^) 2 (^t :i ) 2 < lGeV, where we take 
L = 10 fm and a = 4. In Fig. 7, the time evolution of chiral condensate with the expanding 
geometry is depicted. It is seen that the behavior of damped oscillation is realized quicker 
than the case without spatially expansion. In Figs. 8 and 9, the time evolution of the lowest 
(n = 0) and the first excited (n = 1) pion modes is depicted. They show that the amplitudes 




Fig. 9. The time evolution of the fluctuation Fig. 10. The time evolution of the fluctuation 
mode with n = 1 in the 7r-direction is de- mode with n = in the a-direction is de- 

picted, picted. 
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of these modes are amplified weakly. However, the sigma modes are not amplified even the 
lowest sigma mode (n = 0) as is seen in Fig. 10. 

It should be here noted that the amplification occurs in the late time of chiral phase 
transition both the case of no spatially expansion and of the one-dimensional spatial expan- 
sion. Namely, even when the chiral condensate oscillates around its vacuum value with small 
amplitude, there are amplified solutions of pion modes. The mechanism of amplification is 
clarified in the next section. 

§4. Late time of chiral phase transition 

In this section, we investigate the time evolution of the chiral condensate and quantum 
meson modes based on Eqs.(2-22) and (2-23) or (2-27) and (2-28) in the late time of chiral 
phase transition. Especially, we concern with the amplification of quantum meson modes, 
which was shown in the previous section. It is shown that possible mechanism may be an 
amplification by the forced oscillation as well as the parametric resonance. 

4.1. Linear approximation around static configurations without expansion 

Hereafter, we assume that fluctuation modes of the pion fields are identical one another as 
was assumed in the previous subsection. Also, the condensate depends on time only. When 
the explicit chiral symmetry breaking term, h, is small, the static solutions in Eq.(2-22) and 
(2-23) are given as 

1 3 ~ h 
^ = ^Va M °-2MJ' 

u%t)=u?(t) = -^e- iE Z t , (4-1) 

where E% = \fk 2 + M%. The Fourier mode 1 / a/2~E* has been determined by the normal- 
ization condition (2-14). 

Let us investigate the time- dependent solutions around the above static configurations. 
The condensate <p and quantum meson modes can be expanded around the static solu- 
tions of Eq.(4-1) : 

<Po(t) = y?o + 6<p(t) , 

= + • (4-2) 

Here, we consider the late time of the chiral phase transition. Then, \5(p(t) | is small compared 
with the vacuum value tp . Further, we assume that \5u^\ <C \u^ s \. Also, in general, the 
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fluctuation G given in (2-4) is small compared with </?q. Since G can be written in the form 
in (2-16), we conclude that the following relation should be satisfied : 

1 



(2tt) 3 

£oM*)l» 



d 3 k \v£ 
1 



d*k(u k a (t)5u k a (t) +u k a (ty5u k a (t)) 



(4-3) 



Under the above approximation, the equation of motion for 5ip(t) and Su^(t) are obtained 
from Eqs.(2-22) and (2-23). For the condensate 5<p (t), 



(d t + M*)6<p(t) = 
is obtained. A solution of (4-4) is written as 

5(p(t) = -5acos(M t + (f>) , 



(4-4) 



(4-5) 



where 5a and are constants. We rewrite hereafter a = and 1 ~ 3 into a and it. Let us 
introduce the following dimensionless variables : 

5u*(t') = y/2M~5u*(t) , 7 t' = M a t , 



7 2 (£; 2 + M 2 ) 



Ml 



7 2 (fc 2 + M 2 ) 



h = V ° 



6a h 1 



Ml 



8a 



F„ = 



3 P + y? 

m 2 y /4 5o- 



k 2 + Ml 
Ml 



3M 2 U 2 + M 2 



1/4 



V?0 

5a 

<Po 



(4-6) 



where we put 7 = 2. If we set up <fi = without the loss of generality, then the equations of 
motion for the meson modes derived from (2-23), (4-1) and (4-2) can simply be expressed as 

d 2 



dt' 2 



+ u a [l- h a cos(7t') 



8ul{t') = F a cos( 7 0e w , 



(4-7) 



where a = a or n. If F a (F n ) is negligible, then Eq.(4-7) is reduced to Mathieu's equation. In 
this case, the existence of the unstable solution for 5u^(t') may be expected. This phenomena 
is well known as parametric resonance in classical mechanics. On the other hand, the forced 
oscillation may be realized by the effect of F a (F n ) even if the model parameters do not offer 
the unstable regions for the parametric resonance. 
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4.2. Linear approximation around static solutions with expansion 

In this subsection, we investigate the case with the expanding geometry. We derive the 
equations of motion with expansion using the same formalism in §4.1 in the late time of 
chiral phase transition. The solution of (2-28) is given as 

ftj^W'ME^, (4-8) 



where E^ = y (A; 2 + 1/4) /r 2 + k\ + M 2 . We expand the proper-time dependent variables 
<Po(t) and ut\r) in (2-26) around ip and the above solution (4-8), respectively : 

^o(t) = y?o + S(po(r) , (4-9) 
v^(t)=v^'{t)+Sv^(t). (4-10) 

Corresponding to (4-3), the following approximations are adopted : 

VoSMr) » \- j^y 3 Jfd^Jdk n {u^\r)5u^\r) + u^\t)5u^\t)) . (4-11) 

Under the above approximation, the equations of motion for 5<po(r) and 5u]z\t) are obtained 
from Eqs.(2-27) and (2-28) as 

(d 2 T + U T ) + M CT 2 (r)} 5 Vo (r) = , (4-12) 

d 2 T + S^i^ + k\ + Ml + Ay? • 5Mr?) 6uW(t) = -XMr) • ^(tR 1 ^-^ , 

+ + fcl + Ml + • ^o(r)) 5uM(t) = -^o(r) • ^(t)^-^ . 

(4-13) 

Here, we have omitted dE^/dr because we will replace r into r c in the second term in 
(2-25) in order to stress the effect of friction due to the spatial expansion, where r c may be 
taken as the proper time when the hadronization occurs. Thus, by putting r into r c of the 
second term in the left-hand side in (4-12), the amplitude of chiral condensate around the 
vacuum value reveals a behavior of a damped oscillation. As is well known, this equation 
has three-type solutions. For (1/t c ) 2 < 4M 2 , we obtain 

6(Pq(t) = -5ae' T/2Tc cos(cur + 0) , (4-14) 
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where 5a and are determined by the initial conditions and 



CO 



Mn o • 



(4-15) 



Here, we adopt = for simplicity. We define the following dimensionless variables : 



7 



V±Kr? - 1 ' 



0J n 



7 2 ((k 2 v + l/4)/r c 2 + k\ + Ml) 



Xip 2 



5a 



_ _ ____ S ± F ^tM.J^yM. 5 ! 

~3 (k* + l/4)/T* + ki + M* <p ' w 3^ V ^ ^ ' 



^ (fcg + l/4)/r c 2 + ^ + Ml <p ' 



h w = - ■ 



\ip 2 



(4-16) 



where a = a or ir. By introducing the above dimensionless variables, 5ipo and the equation 
of motion for quantum meson modes in (4-13) are recast into 



5(p (t') = -5ae'^' 008(70 , 

^ + cu 2 a [l - K cos( 7 t / )e- /1 *'] ) ^i k] (0 = F a cos^O^' 



-itJat' 



(4-17) 
(4-18) 



The unstable regions in which the amplified solutions of (4-7) and (4-18) exist will be inves- 
tigated in the next section numerically. 

§5. Unstable regions for quantum meson modes 

In this section, we numerically show the unstable regions in which the absolute values 
of the amplitudes 5u derived from solutions of the equations of motion for quantum meson 



Fig. 11. The unstable regions are depicted in the usual Mathieu equation. 
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modes in (4-7) without spatial expansion and in (4-18) with spatially expansion, respectively, 
are amplified. This phenomena are seen in the lowest and the first excited pion modes as 
was demonstrated in the previous numerical calculations. 

Hereafter, we omit the subscript a. In Fig. 11, the unstable regions, which are represented 
in the dark area on the figure, for the usual Mathieu equation are depicted, which corresponds 
to the case F a = in (4-7). The horizontal axis represents u 2 h/2 and the vertical axis 
represents u 2 , respectively. On the other hand, for F a ^ 0, the unstable regions are added 
to those given in the usual Mathieu equation because of the effect of the right-hand side 
in (4-7) which reveals the effect of the forced oscillation or the beat. Here and hereafter, 
we obtain the unstable regions, in which the amplification of the quantum fluctuation mode 
functions occurs, as follows : We performed the time-integration in a certain time interval 
of the magnitude of the quantum meson mode function, which means a kind of the time 
average. Then, we compare a time-integration in a certain time step with the proceeding 
one. If the value of time-integration is larger than the last one for all time steps under 
consideration, then we decided that this mode is unstable. In Figs. 12 and 13, the cases 
F a = 1.0 and 3.0 are shown, where the time-integration is taken in 10 fm/c in order to judge 
the unstable or stable solutions up to 100 fm/c. From these results, it is obviously found 
that the unstable regions become wider due to the effects of the forced oscillation and/or 
the beat induced by the oscillation of the chiral condensate. 

For the case in Eq.(4T8) where the friction term appears in the equation of motion 
because the system has the spatially expanding geometry, the unstable regions are depicted 
with F a — 0, 1 and 3 in Figs. 14, 15 and 16, respectively. The parameter \x is determined 
by adopting M a = 500 MeV and r c = 10 fm/c . The time interval for the time-integration 
to determine the unstable regions in our approach is taken as rather short time, 4 fm/c, 
because the unstable regions are disappeared by strong damping effect for long time past. 
We calculate the time evolution up to 40 fm/c. In Fig. 14 with F a = 0, the result for the 



F=l 



F=3 



Fig. 12. The unstable regions are depicted for Fig. 13. The unstable regions are depicted for 
the solutions in Eq.(4-7) with F a = 1.0. the solutions in Eq.(4-7) with F a = 3.0. 
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Fig. 14. The unstable regions are depicted in the Mathieu equation with friction term in (4-18), 
in which we set F a = 0. 





F = 1 


H,„ ' .. 







u 2 h/2 uPhjl. 

Fig. 15. The unstable regions are depicted for Fig. 16. The unstable regions are depicted for 
the solutions in Eq.(4-18) with F a = 1.0. the solutions in Eq.(4-18) with F a = 3.0. 

unstable regions is well known because this case corresponds to the case of the Mathieu 
equation with a friction term. In Fig. 15 and 16, we introduce the effect of F a in the 
realistic situation in Eq.(4-18) which gives the forced oscillation or the beat. However, these 
effects only give a small modification. The reason is that the friction works strongly in this 
model. As a result, the amplification of the quantum meson modes with spatially expanding 
geometry is not so strong, although the amplification appears. This results mean that the 
effect of forced oscillation or the beat ceases to work for spatially expanding geometry. 

§6. Summary 

We have demonstrated that the amplitudes of quantum pion modes are amplified even 
in the late time of chiral phase transition by the mechanism of the forced oscillation as well 
as the parametric resonance in the framework of the 0(4) linear sigma model. The basic 
equations of motion have been derived in the time-dependent variational approach with the 
Gaussian wavefunctional, in which the equations of motion of the chiral condensate and 
quantum meson mode functions have been formulated in the form of the Klein-Gordon type 
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equation and the Liouville-von Neumann type equation of motion, respectively, in both the 
cases of the uniform and the one-dimensional spatially expanding systems. These equations 
of motion have been solved numerically, and it has been shown that the amplified solutions 
for the quantum pion modes with low momenta exist really. 

We have investigated the mechanism for this amplification phenomena, and have pointed 
out that there is a possibility of the resonance by forced oscillation or the beat as well as the 
parametric resonance induced by the oscillation of the chiral condensate around its vacuum 
value, while the effect of the forced oscillation ceases to work due to the strong friction 
for the case that has the spatially expanding geometry. Also, by numerical calculation, we 
have concretely given the parameter regions in which the unstable solutions for the quantum 
meson mode functions exist. Of course, it is necessary to know the initial conditions in 
order to judge whether the amplification of quantum meson mode occurs or not, namely, the 
parameters are in the unstable region or not, in the realistic relativistic heavy ion collisions. 
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